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AN ANALYSIS OF LAMINAR FREE-CORVECTIOR FLOW ARD BEAT TRANSFER 
ABORT A FLAT PLATE PARALLEL TO THE DIRECTIOH 
OF THE GENERATING BODY FORCE 
By Simon Ostrach 


SUMMARY 

The free -conyect ion flow and heat transfer (generated hy a tody 
force) about a flat plate parallel to the direction of the body force 
are formally analyzed and the type of flow is found to be dependent on 
the Grashof number alone . For large Grashof numbers (which are of 
interest in aeronautics), the flow is of the "boundary-layer type and the 
problem is reduced in a formal manner, which is analogous to Prandtl' s 
forced -flow boundary -layer theory, to the simultaneous solution of two 
ordinary differential equations subject to the proper boundary conditions . . 

Telocity and temperature distributions for Prandtl numbers of 0.01, 
0.72, 0.733, 1, 2, 10, 100, and 1000 are computed and it is shown that 
velocities and Eusselt numbers of the order of magnitude of those 
encountered in forced-convection flows may be obtained in free -convection 
flows . The theoretical and experimental velocity and temperature distri- 
butions are in good agreement. 

A flow and a heat -transfer parameter, from which the important phys- 
ical quantities such as shear stress and heat -transfer rate can be com- 
puted, are derived as functions of Prandtl number alone. Comparison of 
theoretically computed values of the heat-transfer parameter with values 
obtained from an appro xima te calculation and experiments yielded good 
agreement over a large range of Prandtl number. Agreement between the 
theoretical values and those obtained from a frequently used ssmiempiri- 
cal heat-transfer law was good only in restricted Prandtl number ranges 
(depending on an arbitrary constant) . 


INTRODUCTION 

Two important types of fluid flow problems involving heat transfer 
are those of forced and those of free convection. By forced-convection 
flow is meant flows maintained mechanically as, for example, by a pres- 
sure drop or an agitator. Free -convection flow, on the other hand, 
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results from the action of "body forces on the fluid, that is, forces 
which are proportional to the mass or the density of the fluid. The 
flow is generally produced in the following manner: Consider, for 

example, a fixed object (such as a plate) in a quiescent fluid subject 
to a body force. When the plate is at the same temperature as the sur- 
rounding fluid, the body forces acting on the fluid are in equilibrium 
with the hydrostatic pressure and no flow ensues in the steady state. 

If a temperature gradient normal to the body force is Imposed by heating 
(or cooling) the plate, there will exist a defect (or excess) of body 
force because of the decreased (or increased) density, with the fluid 
closer to the plate having the greater defect (or excess) than that away 
from the plate. This unbalance of the forces causes the fluid to be 
accelerated with the particles nearer the plate moving more rapidly than 
those farther from the plate. Free -convection flow has usually been con- 
sidered to be generated in a gravitational field where the previously 
mentioned defect or excess of body force was the Archemedian (buoyancy) 
force. However, since centrifugal farces are also proportional to the 
fluid density, free -convection flows can also be set up by the action of 
such forces. (See reference 1 for a more explicit discussion of the 
development of free-convection flowB by centrifugal forces.) 

Free -convection flows produced by centrifugal forces are now of 
practical Importance in aeronautics because many aircraft propulsion 
systems contain components (such as gas turbineB and helicopter ram jets) 
which rotate at high speeds and in which heat is being transferred. The 
method of free-convection cooling of gas-turbine rotor blades where the 
centrifugal forces create a free-convection flow of the coolant in the 
blade passages is an example of a practical application of the free- 
convection phenomenon in aeronautics . Also, free-convection flow due 
to centrifugal force is superimposed on the flow through helicopter ram 
jets and on the flow of cooling air in hollow rotor blades of air-cooled 
turbines and, under proper conditions, can appreciably Influence the 
resultant flow and heat transfer. 

As a simplification of the many free-convection problems which are 
now of some consequence in aeronautics, consideration is here given to 
the special case of free-convection flow about a flat plate parallel to 
the direction of the generating body force. The experimental and theo- 
retical considerations of Schmidt and Beckmann (reference 2) concerning 
the free-convection flow of air subject to the gravitational force about 
a vertical flat plate constitute the most complete treatment of this 
subject up to the present time. Eckert (reference l) as well as others 
has further verified and extended the experimental results of Schmidt 
and Beckmann, and Schuh (reference 3) has extended the numerical calcu- 
lations by computing the velocity and temperature distributions for 
several Prandtl numbers different from that for air. However, all the 
'theoretical work in these references is based on the incompressible 
equations in which the density (or temperature) variation is introduced 



2370 


KAjCA TKf 2635 


3 


in tlie "buoyancy term, alone. Various terms are omi tted from the equations 
at the start on the basis of either intuitive arguments or no arguments 
at all. Although a theoretical development made in such a manner led to 
good final results, the significance of all the important factors assoc- 
iated -with the free -convection flow phenomenon is not obtained from such 
an analysis. 

The problem of free -convection flow as produced by a body force 
about a flat plate in the direction of the body force was studied at 
the BACA Lewis laboratory and is treated in a fo rmal and more general 
manner herein. The method used is somewhat similar to that used in 
reference 4 wherein consideration was given to the free -convection flow 
at high Grashof numbers in a horizontal cylinder which had a variable 
surface-temperature distribution. The application of this method to 
the present problem leads to a development which is analogous to Prandtl's 
treatment of high Reynolds number forced -convection flows. Although the 
final equations obtained by this method are the same as those of Schmidt 
and Beckmann, this more general approach not only clearly demonstrates 
the significance of all the important parameters and assumptions and 
hence leads to a better understanding of this type of flow but also 
indicates the quantitative limitations of the theory. In addition, the 
numerical solutions of references 2 and 3 are herein extended to cover 
a more complete range of variables . The new calculations yield infor- 
mation on the free -convection flow for Prandtl numbers corresponding to 
those of liquid metals, gases, liquids, and very viscous fluids. 


A&ALXSIS 

Statement of Problem and Basic Equations < 

The steady -state equations expressing the conservation of mass, 
momentum, and energy for a compressible, viscous, and heat-conducting 
fluid subject to a body force together with an equation of state govern 
the flow and associated temperature distribution about the plate. These 
equations in Cartesian tensor notation are (see reference 5), respec- 
tively, 

(PUJ = 0 (1) 


Su j\ _ j3P 
oxJ dXp 


ax. 


au. 


P TT 


i ax, 


= Pf 


1 + ax. 
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p = p(p,t) 


( 3 ) 

( 4 ) 


(A complete list of the symbols used herein is given in appendix A.) 

For the two-dimensional case, equations (l) to (4) represent a system 
of five equations in the five dependent variables TJ^, TJg, P, P, and T. 
For later use, equation .(4) can be written 

dp = p(E dP - p cLT) (4a) 


where K and p are the coefficients of isothermal -compressibility and 
volumetric expansion, respectively (see reference 6) . In addition to a 
general state equation, such as is given in equations (4) or (4a), it 
will be convenient at times in the discussion to refer to some specific 
state equation. To this end, the equation of state for an ideal gas 

P = pET (4b) 


will be used. 

Particular consideration is here given to the two-dimensional free- 
convection flow about a semi-infinite vertical flat plate . The X-^-axis 
of the coordinate system is taken along the plate and the Xg-axis , nor- 
mal to it. Ho distinction is made as to the specific type of body force 
acting, f<pr example, gravitational or centrifugal, but the force is 
assumed to be acting in the vertical direction only (that is, parallel 
to the plate) . Centrifugal and Coriolis forces which are connected with 
flows on curved paths and with rotating systems generally vary with 
position and velocity. However, in order not to mak e the analysis 
unduly complicated, the body force is- taken to be constant. 

Tn order to define the problem clearly, a choice must still be made 
of the position cf the origin of the coordinate system. Before making a 
definite decision on this point, note that for constant plate tempera- 
tures there are four permutations of the body -force direction (either 
upward or downward) and the plate thermal condition (either heated or 
cooled) which will lead to free -convection flows. Once the position of 
the edge of the plate, which is also to be the origin of the coordinate 
system, is decided, there are two combinations of the body-force direc- 
tion and plate thermal condition that will yield flows which proceed 
away from the edge . It is this type of flow that is amenable to the type 
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of analysis to t>e made here. Uhls point will "be more fully discussed 
subsequently. If the edge of the plate (recall that a semi -infinite 
plate has hut one edge) is taken at the bottom of the plate (that is, 
the plate extends to +<*> in the -direct ion), the two combinations lead- 
ing to flows in the proper direction (upward in this case) are, respec- 
tively, the body force acting downward with a heated plate and the body 
force acting upward with a cooled plate. The equations developed for one 
of the cases reduce directly to those for the other. The remaining two 
permutations, namely, the body force acting downward with a cooled plate 
and the body force acting upward with a heated plate, would yield flows 
which proceed downward or toward the edge of the plate if this edge were 
taken at the bottom of the plate. This type of flow would violate a 
physical condition of the problem which states that the flow starts at 
the plate edge. The latter combinations hence will not be considered 
further. 

Because the two acceptable configurations can be reduced essentially 
to one, for the development' to be given here, the origin of the coordinate 
system -trill be taken- at the bottom of a heated plate, with the body force 
acting downward. The assumption is now made that the viscosity and heat- 
conductivity coefficients are functions of the temperature alone and obey 
the following laws : 



The choice of the body -force direction together with equations (5) alters 
equations (2) and (3) so that they become 



( 7 ) 


Note that the only nonzero component of the body force is the -component . 

Boundary Conditions 

The boundary conditions associated with the given problem are that: 



6 


HA.CA IN 2635 


(a) The fluid must adhere to the plate (the no-slip condition of 
viscous flows) and the plate must he a streamline, or mathematically, 

XTid^O) = U 2 (I 1# 0) = 0 (8) 

(h) The temperature of the fluid at the plate must he equal to the 
plate temperature, that is, 

T(*1,0) = T 0 (9) 

(c) The velocity at large distances from the plate must he 

undisturbed, or 

TTlUi,”) = 0 (10) 

(d) Th e temperature at large distances from the plate must he equal 
to the undisturbed fluid temperature, or 

TfX-L,-) = T. _ (11) 


Simplification of Equations 

Let a small quantity e now he defined as 

e=P(T 0 -T.) (12) 

•which is a measure of the magnitude of temperature variation in the flow- 
field. The coefficient of volumetric expansion p is. generally of the 
order of magnitude between 10“^ and 10 ^ (see table 15 of reference 7, 
for example) and for gases, p = l/T. (Thus, for gases, if p is taken 
to he constant, e = (T 0 -T«,) /T» ; that is, a is the relative temperature 

difference.) The coefficient p will he assumed constant. Because in 
the steady state flow ensues only when there is a temperature variation 
in the fluid, the free -convection velocity should then depend directly 
on e, e-nrl the variations in pressure and density (from the static, 
e = 0, case) due to the temperature differences should also depend on e 
Thus . 


& 

II 

CO 

\ / 

(13) 


p = p s 

+ P 006 CT 

( 11 ) 


P = Ps 

+ PcoSCp 

(15) 

L* 

T = T»(l + 60) 

(16) 
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where “ f X denotes the X]_ -component of the "body force per unit mass, 
u-^, cr, <p, and 6 denote dimensionless functions (which, in general, can 
he functions of e ) , 2 is some characteristic length (for example, the 
distance from the edge of the plate to the point of interest), P g and 
P s are the pressure and the density, respectively, for the static case, 
(UjL = 0 or e = 0), and P* and denote constant values of the 

pressure and the density (that is, the values if no force field were pre- 
sent) defined hy the state equation (in the case of a gas, in particular, 
P<» = PoBTco) . Because there is no characteristic velocity associated with 

the type of flow under consideration, the velocity is dimens ionalized hy 
the factor given in parentheses on the right side of equat ion (15). This 
spec ific factor is chosen because it leads to si mpler equations. ’ 

In order to determine the static quantities, it will at first he 
convenient to consider the particular case of a gas. The problem is 
then considered with the temperature uniform throughout the flow field 
at the value (therefore there will he no flow and = 0). For 

this situation, equations (4b) and (6) become 

P Q = p s ET„ (17) 


and 


BP 

a 

BXi 




+ P s f X = 0 


Bp 

e 

BXc 


> 


( 18 ) 


= 0 




(it should he noted that equation (18) expresses the physical fact pre- 
viously stated that the body force and hydrostatic pressure are in equil- 
ibrium for the static case.) Substitution of equation (17) into equa- 
tion (18) leads to 


p a = So exp - 


-2L 

BT. 


Xi 


(19) 


and equation (19) together with equation (17) and the equation defining 
Po» and p* yields 


fV P. 


P s = ftn exp 1 - XJ = P» exp I- -p 


X 




( 20 ) 


If the exponential in equation (20) is expressed in terms of its series 
expansion, that equation becomes • 
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p *? X T 

■?r Zi 


+ 



( 21 ) 


A computation of the second term, in the brackets of equation (21) for the 
case of air under normal conditions -with f j = g and the fact that 
is of the order of magnitude 2 show that p^z/P^ ~ 10“^l/f oot . For 
the type of problem under consideration, 2 "will always he of unit order 
of magnitude so that eve*n if the "body force f X represents a centrifugal 
force many times that of gravity, the inequality pu m « 1 may 

still he satisfied. Thus, in the subsequent development it will he 
assumed that p s s p^,. This assumption which was Justified hy the com- 
putation for the case of a gas is expected to he reasonable for other 
fluids as well. The physical interpretation of this assumption is that 
under static conditions (e = 0), the density (or pressure) is not 
affected hy the force field. 


In order that all quantities in the following equations he dimen- 
sionless, it is further necessary to define x^ = X-j/2, where the x^ 
are now dimensionless space coordinates. Substituting these new coor- 
dinates along with equations (13) to (16) into equations (l), (6), (7), 
and (4a) and noting equations (18) and that p B £ yield, on neglec- 

tion of terms of higher order in e compared with those of order e, 


c>U]_ Sug 
^X]_ Sxg 


(22) 


Gr |ut ^ + u 2 = Au x - NGr ^2- - cp 


du- 

Su 


ch:-. 


Sue 


Gr l u i &r + ^ ~i = AU2 " ** 


bo 
bx o 


(23) 

(24) 


„ _ ^ f be be 

Gr Pr [ u*i - — + Uo - — 
\/ Sx x ^ dx 2 


^ GrPr [ U ^ +U 2 bt) + A6 (25) 


d<p = E ^dcr - PTa, d0 


(26) 


where P»/fiofj2 = UGr and the Grashof number Gr and the Prandtl number 
Pr are defined as 


(hr = 


pJ^X^e. 


and 


Pr = 




£ 

to 

CV] 
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Physically speaking, the Grashof number represents the ratio of the tody 
forces to the viscous forces. The free-convection flows of interest here 
are those associated with large Grashof numbers . The factors K and p 
in equation (26) may well he taken to he constapts (see reference 6). 

The boundary conditions (equations (8) to (ll)) in terms of the new 
dimensionless variables are 

ui(xi,0) = u^XljO) = 0 (27) 

e(x l' 0) = pfe ( 28 ) 


Ul(zi,») = 0 ' (29) 

0(x 1; «) = 0 (30) 

Thus, to a first approximation, equations (22) to (26) together with the 
boundary conditions replace the original equations and boundary conditions 
(Note that for gases, pT*, = 1.) 

The prime assumption that has been made in this analysis is that the 
higher-order terms in e were negligible, which implies that e is 
small and, consequently, that the temperature difference or p is 
moderately small. It is a consequence of this assumption alone that the 
basic equations were simplified to equations (22) to (26) wherein the 
viscosity term in the energy equation is neglected and the only coupling 
of the momentum and energy equations occurs by means of the body-force 
term in equation (23) . As a result of this assumption, the variations 
of the viscosity and heat -conductivity coefficients with temperature are 
also negligible. Without any discission, the authors of reference 2 
start directly from simplified equations of the same form wherein the 
pressure terms in the energy equation were also neglected. In refer- 
ence 3 some intuitive arguments are given to justify the simplified 
equations . 


It is now convenient to revert to the more familiar notation where 
x = Xp, y = x 2 , u = up, and v = U 2 » Equation (22) implies the existence 
of a stream function Uf such that 


u 


Sy 7 


and 


v 


5y 

dx 


- Y x (31) 


where subscripts denote differentiation. Applying equation (3l) to 
equations (23) to (25), yields, respectively, 
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Gr ^7 ~ ^ Y y*zy " ^x^yy + 

_ - & (£*x) = -Y /22 -f ¥ 3*37 + ^y (33) 

Gr Pr = Y y e x “ ^x?y “ 5^1 ( ¥ y°x “ Y x a y) ( 34: ) 

The boundary conditions (equations (27) to (30)) 1)800100 

¥ y (x,0) =¥ x (x,0) = 0 (35) 

e(*,°) = ^ (36) 

¥y(x,~) = 0 (x,») = 0 (37) 


Equations (32) to (34) and equation (26) form, the system of equations 
for the four unknown functions Y, 9 , <p, and a of the problem. The 
system is nonlinear, and therefore further simplification of the equations 
■would he desirable. Just as in the case of forced-convection flowB where 
the Eeynolds number determines the type of flow or, in math em atical terms, 
the type of solution, the Grashof number is the prime factor for free- 
convection flows. For the case of small Grashof number, it can be seen 
from equations (32) to (34) that a perturbation in the small parameter 
Gr ■will yield a system of linear equations . For Grashof numbers of unit 
order of magnitude, no further important simplification can be made and 
the solutions would have to be obtained numerically. For the other lim- 
iting case, that of large Grashof numbers (which is the case under con- 
sideration herein), it would, at first thought, appear that some simpli- 
fication could be obtained by performing a perturbation in the small 
parameter 1 /Gr. However, this would then imply that the term containing 
the highest -order derivatives (the left term in equations (32) to (34)) 
could, among others, be neglected. (This argument would also imply that 
the body -force term <p in equation (32), which is essentially causing 
the flow, could also be neglected.) The omission of the highest -order 
derivatives from consideration, however, would lead to solutions which 
would not satisfy all the boundary c.onditions . Problems of this type 
are referred to as singular perturbation problems. For further dis- 
cussions of singular perturbation problems see references 8 and 9. 

Equations in which a small parameter multiplies the highest-order 
terms are said to be of the boundary-layer type, because in order for 
solutions which satisfy all the boundary conditions to be obtained, 
the highest-order terns must be considered near the boundary. This 
fact implies the existence of a thin region, called the boundary layer, 
wherein the functions vary rapidly from the value at the boundary to 
that in the flow outside this layer. The conclusion to be drawn from 
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the preceding discussion is that for large Grashof numbers the flow 
is of the hound ary-layer type. Schmidt and Beckmann (reference 2} 
also made the boundary -layer assumptions in their theoretical develop- 
ment, and these assumptions were justified on the basis of their experi- 
mental observations. The Grashof numbers for their experiments were of 
the order of 8X10 6 . 

In view of the fact, previously discussed, that highest-order deri- 
m vatives of each dependent variable as well as of those terms of physical 

-o importance (as, for example, the body -force term) must be retained in the 

boundary layer, it is convenient to make both sides of each of the equa- 
tions of the same order in Gr. In this way, as will be shown, the 
equations will be simplified further. It is thus convenient to make the 
following transformations in the system of equations (32) to (34) and 
(26) and then to retain only the dominant parts (that is, those multi - 
plied by Gr to the highest power) of each individual term. 

Let y = Gr 1 " y, Y = Gr 3 f, cr = Gr* 3 a, <p = qT , and 9 = 6. Then 
equations (32) to (34) and (26) became 

Gr s+3r-ly___ _ Gr -1 ~ c Gr 2s+2r ( y_Y^ - ?J__) + TT GrV (38) 

_ Gr s+2r- 1 y = Gr 2s+r (.fLf^ + + ff Gr t+x Z f (39) 


I 2r_1 e n - GlJ3+r <Vx - W - Grtw+r (T r?. - AV ^ 

dtp = KP 00 Gr t do - pi^dS • (41) 

It can now be seen that by proper choice of r, s, and t a transforma- 
tion of the type given provides a means for making the important terms 
in the differential equations of the same order in Gr. Thus if r = l/4, 
s = -3/4, and t es -1, equations (38) to (41) become 

Y yyy ^ = Y y Y xy “ Y 3?yy + N a x 
_l 

WoL =0(Gr 2 ) £ 0 (43) 

y 

®yy = ** % ? X " *X®y>' 
dtp + pG^d0 «= 0 


( 45 ) 
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More generally, if If Is very much different from, unit order of magni- 
tude, a value of t can always "be chosen (depending on U) such that 
equations (42) to (45) are obtained. (For any negative t less than -1, 
the last term, cf equation (42) will also disappear.) 

There are now several important points to he .discussed concerning 
the transformation Just made and the resulting simplified equations. 
First, it should he noted that the transformation is merely a fo rmal 
expression of the houndary-layer assumptions first made hy Prandtl and 
hence the solutions will he asymptotic for large G-r. Second, the second 
equation of motion here also reduces to state that the pressure across 
the "boundary layer is constant. Third, the pressure terms in the energy 
and state equations are here found to he negligible. This fact verifies 
a priori assumptions made hy others from the physics of the problem. 
Finally, note that integration of the general state equation (independent 
of pressure) as now given hy equation (45) leads to 

<P + 0TJF = 0 (46) 

where the constant of integration has been taken as zero without any loss 
of generality. For the particular case cf a gas, p = l/T^ bo that equa- 
tion (46) becomes 

<p + 6 = 0 

The boundary conditions (equations (35) to (37)) now can he written 

Y f (x,0) = Y x (x,0) = 0 (47) 

°) = (48) 

Yy(x,“) = = 0 (49) 

If now it is assumed that ct z = 0 in equation (42) since consideration 

is here being given to a flat plate, and If cp is eliminated frcm equa- 
tion (42) hy use of equation (46), there results the system of equations 

+ - Vjy < 50 > 

0yy = ^(^y^x “ ^X®y) (51) 

Thus the problem has been reduced to the solution of the two simultaneous 
partial differential equations (equations (50) and (51)) subject to the 
boundary conditions (equations (47) to (49)). 
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Final simplification of the equations is made by application of the 
so-called similarity transformation of boundary-layer theory. Thus, let 

^ x (52) 

(4x) 4 

and 

3 

y = (42:) 4 F(n) (53) 


0 E= 


5M 

(3T ro 


(54) 


Then equations (50) and (51) are reduced to the following ordinary dif- 
ferential equations: 

F" ' + 3FF' « - 2F ,S +^H = 0 (55) 

H* 1 + 2PrpF H‘ = 0 (56) 

•where the primes denote differentiations with respect to T| . The reci- 
procal one -fourth power similarity as given in equation (52) is charac- 
teristic of free -convection flows just as the reciprocal square-root 
type is characteristic of the forced-convection flows. The boundary 
conditions become 


F'(0) = F(0) =0 (57) 

H(0) = 1 (58) 

F*(») = H(-) = 0 (59) 

The use of a transformation like equation (52) essentially specifies an 
additional boundary condition, namely, that the conditions to be satis- 
fied at y = ® (or T} *= «°) should also be satisfied at x = 0. , It is 
for this reason that the flows previously discussed which would flow 
toward the edge (downward) are not amenable to this type of analysis, 
for such flows would violate this additional condition, which essentially 
states that the boundary-layer development starts at the edge of the 
plate. 


Solution of the Boundary -Value Problem 

The solutions of the simplified equations (55) and (56), satisfying 
the boundary conditions as given by equations (57) to (59), were obtained 
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by use of an IBM Card Programmed Electronic Calculator. A detailed 
account of the procedure followed in the determination of the unknown 
functions is presented in appendix B by Dr. Lynn IJ. Alters. The 
functions P and H together with their derivatives are given in 
table I for Prandtl numbers of 0.01, 0.72, 0.733, 1, 2, 10, 100, and 
1000. Even though the Prandtl number for air is taken as 0.72 in this 
report, the solutions for Pr = 0.733 were also computed and are presented 
as a check with the Schmidt-Beckmann calculations wherein the value of 
Prandtl number of 0.733 was used. The particular values of the Prandtl 
numbers given were chosen to correspond to those for liquid metals, gaseB, 
liquids (such as water and oil), and very viscous liquids (such as glycerin 
or oils at very low temperatures). 


EESUITS 

Telocity and Temperature Distributions 

By means of the various transformations made in the analysis it can 
easily be verified that 


2 Vp(T 0 -T o )f 2; I 2 V^ 


= I“(ti) 


'T-OLN 


where 


f z (T 0 -^)f3 ' 


Equations (60) to (62) relate the physical quantities to the dimension- 
less functions P and H which are now known. The dimensionless velo- 
city and temperature distributions as given by equations (60) and (6l) 
are presented in figures 1 and 2, respectively, as functions of T) for 
the various values of Prandtl number. The computations made here agree 
with those for Pr = 0.733 as given in reference 2 up to the third sig- 
nificant figure. Por Pr = 10, 100, and 1000, the present results agree 
in general with those of reference 3. Since only curves are presented, 
in reference 3, the precision of the agreement cannot be stated. 
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OZhe maxim u m values of the dimensionless velocity distributions occur 
at larger values of the argument tj as the Prandtl number decreases and 
the velocities decrease with increasing Pr. It should also he noted 
that the dynamic and thermal boundary -layer thicknesses can be estimated 
from the abscissas of figures 1 and 2, respectively, and that for Pr» 1 
the velocity boundary layer is much thicker than the thermal boundary 
layer . 

Phe occurrence of fj (or Gr^- as given by equation (60)), which 
may be very large for flows generated by centrifugal forces, in the 
denominator of the ordinate implies that velocities of appreciable mag- 
nitude can be associated with such free -convection flows. In particular, 
if f j = 10® feet per second squared, which is a reasonable conservative 
figure for present-day rotating systems, e= 0.2 {which is within the 
limits of the theory presented herein), and arbitrarily X = 0.25 foot, 
then the maximum velocity attained at a Prandtl number of 0.72 is approx- 
imately 125 feet per second. Phis value of the maximum velocity could, 
of course, be doubled or even tripled under the proper conditions. One 
limitation to a calculation of this sort, as can be seen by comparison 
of the denominators of the left and middle terms of equation (60), 
should be kept in mind; namely, the limiting Grashof number for laminar 
flows . In lieu of a complete stability analysis on this type of flow, 
this limiting value is taken to be 10®, as indicated in reference 10. 
Consideration of this limitation then implies (see equation (60)) that 
for large laminar velocities either must be large or X must be 

small. 


Comparison with Experiments 

Careful experiments of free -convection flows (as generated by gravi- 
tational forces) about vertical flat plates were made by Schmidt and 
Beckmann (reference 2) in which velocity measurements at various points' 
along the plate were made by means of a quartz -filament anemometer and 
the temperature measurements were obtained by means of manganese- 
constantan thermocouples. Eckert (reference l) performed similar experi- 
ments in which the measurements were made by means of a Zehnder-Mach 
interferometer. Phe results of both sets of experiments are in good 
agreement with each other, but since the data presented in reference 2 
by S chmi dt and Beckmann appear in more detail, these data will be used 
for comparison with the theory. 

Phe experiments of reference 2 were performed on two different (in 
that the edges were smoothed either symmetrically or not) 12- by 25-centi 
meter plates and on one 50- by 50-centimeter plate. It should here be 
pointed out that the results for the two smaller plates were almost iden- 
tical and that the flow was entirely laminar except near the outer edge 
of the boundary layer where the slight turbulence of the room air dis- 
turbed the measur eme nts somewhat. (Phis effect was also observed by 
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Eckert.) large periodic oscillations of tie flow near the downstream 
edge of the larger plate were observed in addition to the slight turbu- 
lence near the outer edge of the boundary layer. Hence the data from 
the larger plate should not be expected to yield completely satisfactory 
agreement with the laminar theory as presented here. 

Since the physical quantities can be expressed in terms of a single 
variable as in equations (60) and (61), it is to be expected that the 
data taken at the various points along the plates should all lie on a 
single line if the data are correlated according to equations (60) and 
(61). Thus for the smaller plates where (Tq-T») = 95.22° E and 
T„ = 518.68° B, equations (60) to (62) became 


U “F'(tl) 

4 .862-y/x 

(63) 

T - 518.68 > 

95.22 - H( "> 

(64) 

tj = 88.26 -2j— 

(65) 


X 4 

The velocity and temperature distributions are so plotted in figures 3 and 4, 
respectively, as are the curves computed theoretically for Er = 0.72. It 
can be seen that the agreement is in general very good for small values of 
tj and somewhat less satisfactory though still rather good for the larger 
values of tj . The scatter in the range of the larger values of T) is 
believed to be caused by the previously discussed room turbulence. It 
should also be noted that the points farthest away from the theoretical 
are those measured near the leading edge. These points should not, of 
course, be expected to agree too well with the theory since the boundary- 
layer assumptions made in the theoretical development imply that the dis- 
tance along the plate is large as compared with the boundary-layer thick- 
ness . Hence, this assumption is invalid near the leading edge. Schmidt 
anrl Beckmann obtained closer agreement between the theory and the experi- 
ments for the temperature data and poorer agreement for the velocity data 
by basing the kinematic viscosity coefficient in equation (62) on the 
plate temperature rather than on the undisturbed stream temperature as was 
done here . 

For the larger plate, (Tq-TJ = 83.7° E and T ra = 527.14° E so 
that equations (60) to (62) became 

— 2 — = F'(n) 

4 . 522 -^ 


( 66 ) 
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T - 527.14 
83.7 


= H(rj) 


(67) 


Tj *= 83.93 —X— (68) 


The velocity and temperature distributions for this experiment are plotted 
in figures 5 and 6, respectively, and again the theoretical curves for 
Pr s 0.72 are included. In figure 5 it can he seen that for large t] 
the agreement is rather poor, particularly for the data for both s mal l 
and large values of X. The poor agreement for small values of X is 
again due to the theory limitation near the edge of the plate and for 
the large values of X, to the fact that the flow was becoming turbulent 
there . 


Flow and Heat -Transfer Parameters 

In addition to the velocity and temperature distributions, it is 
often desirable to compute other physically Important quantities (such 
as shear stress, drag, heat -transfer rate, and heat -transfer coefficient) 
associated with the free-convection flow. To this end, two parameters, 
a flow parameter and a heat -transfer parameter, are derived in appen- 
dixes C and D, respectively. 

The flow parameter 


1 — F' 1 (0) 

( 4 Gr x 3 ) 4 ( v - ^o/* 2 ) 

is presented as a function of Prandtl number in figure 7. Thus the 
various flow quantities for a given set of conditions can easily be 
computed by application of figure 7 . 

The local heat -transfer parameter 


ETu 


(Gr z /4) 


1 

4 


-H'(0) 


as determined here is given as a function of Prandtl number in figure 8. 

A calculation of the local Husselt number from this equation for 

Pr = 0.72 and Gr^- = 10® yields a value of 63.5, which indicates that 

large heat-transfer coefficients can also be obtained -with free-convection 
f lows . 
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On the "basis of a simplified theory (-that is, "by use of integrated 
momentum and energy equations and assumed velocity and temperature dis- 
tributions), Eckert (see p. 162 of reference 7) obtained the appro xima te 
relation 

1 

Eu 0.718(Pr) 2 

1 “ -I 

(Grj/4) 4 (0.952 + Pr) 4 


The curve representing this equation is also presented' in figure 8, and 
it closely approximates (to within about 10 percent) the curve determined 
by the more exact considerations of this report over the entire Prandtl 
number range. A semi empirical equation as given in reference 11 relating 
the average (over the length X) ITusselt number to the Prandtl and 
Grashof n umb ers -which bas been used in the heat-transfer calculations up 


to the present is 


lTu aT = 0.548 


1 



The constant 0.548 pertains specifically to air; for oil it should be 
0.555 (see reference 12) and for mercury, approximately 0.33 (see ref- 
erence 13) . In order to obtain local values from the average ones given 
by the last equation, it is merely necessary to multiply the average 
values by 0.75. (The determination of this reduction factor of 0.75 is 
discussed in appendix D.) Thus in terms of the local quantities the 
s emi empirical relation becomes 

Eu = 0.411 f(Pr)GrJ 4 


%■■■■- = 0.581 (Pr) 4 

(Gr^/4) 4 

The curve given by this equation is also presented in figure 8 and the 
agreement with the theoretical curve is very good for Prandtl numbers 
near unity, not so good for large Prandtl numbers, and very poor for the 
htohII Prandtl numbers ■ Of course, changes of the constants in the semi- 
empirical relation as previously discussed for the large or small Prandtl 
number cases (oil and mercury, respectively, for example) would cause the 
semiempirical curve to approximate the theoretical curve more closely. 

The values of the heat -transfer parameter obtained experimentally for 
mercury (Pr = 0.03), air (Pr — 0.72), water (Pr =7), and oil 
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(Pr = 75.5, 115, 190, 224, 275, 318, 368, and 442) are here reduced "by 
the factor 0.75 from the average values reported. The value for mercury 
is an average taken of four readings from a curve since this exper ime nt 
■was the only one not reported in tabular form, item figure 8 it can be 
seen that all of the experimental values except those for the oil experi- 
ments are in very good agreement -with the theoretically computed values . 
The data from the oil experiments, though not so good, show reasonable 
agreement (maximum error of approximately 20 percent) -with the theoreti- 
cal curve and good agreement, as is to be expected, -with the semi- 
empirical curve. The difference between the theoretical values and the 
oil experiment results can possibly be due to the fact that the viscosity 
changes in oil are large even for small temperature differences or due 
to the end effects in the measurements . 


CONCLUSIONS ' 

An analysis was made of the free -convection flow about a flat plate 
'oriented in a direction parallel to that of the generating body force 
under the prime assumption that the relative temperature difference is 
small. It was found that the Grashof number was the principal factor 
determining the type of flow and that for large Grashof numbers the flow 
was of the boundary-layer type. The theoretical development was then 
continued to consider only the cases of large Grashof number because 
these are of most importance In aeronautics . - 

Telocity and temperature profiles for Prandtl numbers of 0.01, 0.72, 
0.733, 1, 2, 10, 100, and 1000 were computed on the basis of a constant 
body force and plate temperature and agreement with experiments where 
the fluid was air (Prandtl number of 0.72) was good. It was also demon- 
strated that, velocities and Nusselt numbers of the order of magnitude of 
those obtained in forced-convection oould be obtained In free-convection 
flows . 

A flow parameter and a heat-transfer parameter which are functions 
of the Pr and tl n umb er alone were derived. Calculations of the important 
physical quantities such as shear stress, heat -transfer rate, and the 
like can be computed from these parameters . Talues of the heat-transfer 
parameter obtained from an approximate theoretical development and from 
experiments compared with values computed from the present development 
showed good agreement over a. wide range of Prandtl number (0.01 to 1000). 
It is shown that the commonly used semiempirical relation for the heat- 
transfer coefficient will yield good results only In restricted Prandtl 
number ranges. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 3, 1951 
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SYMBOLS 


The following notation is used in this report: 


(»), L..( n ) 


C' 


X x 

Gr 


coefficients in numerical differentiation and integra- 
tion' formulas 

specific heat at constant pressure 
dimensionless Telocity function 

components of "body force per unit mass, i = 1, 2, 3 
ne gat Its of X-camponent of "body force per unit mass 


Grashof number. 


~ 2* ,3 

poo fjt £ 


Gr x 

8 

H 

h 

K 

k 

l 

m, n 

Iff 

Uu 

Hu 
. av 

P 

Pr 


au/p) 


L ap J 


Grashof number "based on X 

gravitational force per unit t mass (or acceleration 
die to gravity) 

dimensionless temperature function 
"h flat -tr ans fer coefficient 

isothermal compressibility coefficient, -p 

h.eat -conductivity coefficient 
characteristic length 
arbitrary exponents 

a number, defined following equation (26) 

Uusselt number, hX/k 
aTerage Husselt number 
pressure 
Prandtl number 


T 
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E 

r, b, t 
T 


*i 

u 

T 


X 


1 


X 

X 


1 


y 

P 


7 

A 

e 

n 

e 

X 

© 

M- 

v 

P 

a 

T 

<P 

Y 


gas constant 
arbitrary exponents 
absolute temperature 
Telocity components, 1 = 1 , 2 , 3 
dimensionless Yelooity components, i = 1^ 2 , 3 
dimensionless Telocity component in x-direction 
dimensionless Telocity component in y-direction 
Cartesian coordinates, i = 1 , 2 , 3 
dimensionless Cartesian coordinates, i = 1 , 2 , 3 
dimensionless Cartesian coordinate 
dimensionless Cartesian coordinate 


coefficient of Yolumetric expansion. 



ratio of specific heats 
Lap lac i an operator 

relatiTe temperature difference, 0(Tq-T«>) 

similarity -variable 

dimensionless temperature function 

step size used in numerical calculations 

absolute -viscosity 

kinematic -viscosity 

density 

dimensionless pressure function 
shear stress 

dimensionless density function 
stream function 
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Subscripts : 

s 
0 

0* 

Subscript notation 
Superscripts : 

Primes denote ordinary differentiation. 

Bars (as a or y ) denote transformed dimensionless quantities. 


Cartesian tensor and summation subscripts 
denotes evaluation at static conditions (e = 0) 
denotes evaluation at plate surface 
denotes evaluation at undisturbed conditions 
is used to denote partial differentiation. 
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APPENDIX B 

U0MEEICA1 S0L0TI01T OF SIMPLIFIED BOUNDARY -VALUE PEQBLEM 

By Lynn TJ. Alters 


The method Is presented herein by which solutions to the 'boundary- 
value problem 


F' " + 3FF’ ’ - 2F’ 2 + H = 0 • 

(Bl) 

H' * + 3 Pr FH 1 =0 

(B2) 

F(0) = F’(0) = 0 H(0) = 1 


F'(«) = H(») = 0 



were obtained for the cases of Pr equal to 0.01, 0.72, 22/30, 1, 2, 10, 
100, and 1000. This discussion frill enable the results to be clearly 
evaluated and frill perhaps serve as a guide in the numerical solution of 
similar problems. 

Each of the cases of the problem has a solution for a particular 
set of values for F''(0) and H’(0), hereinafter called eigenvalues. 

The basic approach to the problem was to estimate the eigenvalues and to 
integrate out from zero, obtaining functions which satisfied equations (Bl) 
and (B2) at each step. The integration was continued until the functions 
F' and H behaved in a fashion inconsistent with the boundary values at 
infinity; for example, when they became negative or diverged to infinity. 
Improved estimates of the eigenvalues were then made on the basis of the 
results of preceding runs and the process was repeated successively until 
a solution was obtained. 

Modifications Required to overcame specif io obstacles will be dis- 
cussed after sufficient details of the basic procedure have been given. 
Then an evaluation of the accuracy of the numerical results will be made. 

The integration process consists of two parts, a starting phase and 
an extension phase. The starting phase begins with an estimate of the 
eigenvalues F 1 1 (0) and S' (0) and a decision on the step size x to 
be used. It continues with an iterative process of alternately com- 
puting F' 1 1 and H' ' at the first four points and integrating them 
by five-point formulas. This process and that in the extension phase 
are so arranged that the differential equations are satisfied at each 
integral multiple of the step size. 
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The extension phase used preceding data to integrate step by step 
beyond the fourth point. Diagrams of both phases -will be given after a 
few prelim i nary explanations . 


All integration formulas used are based on the same idea. If a 
function, for example, F' ”, is known at five points, there is a unique 
fourth-degree polynomial which agrees with it at these five points . More- 
over, if the successive antiderivatives (integrals) F' ' , F‘ , and F of 
F' ” are known at one point, there are unique fifth-, sixth-, and 
seventh -degree polynomials which are successive antiderivatives of this 
fourth-degree polynomial and which agree with F 1 1 , F' , and F, respec- 
tively, at the one point. It is then a simple algebra problem to deduce 
from the values of F 1 1 ’ at five points and F, F 1 , and F' ' at a 
single point the values of any of these four polynomials at any point. 
These results will approximate the functions F, F 1 , F' ' , and F* 1 ' to 
a degree dependent on step size, the relative positions of the points in 
question, and the magnitude of the fifth derivative of F' " in the 
neighborhood of these points . 

The preceding algebra problem can be presolved in all situations that 
arise in the starting and extension phases of the present problem and 
specific integration formulas may be deduced. These formulas are dis- 
cussed in the next paragraph. 


Let F” 1 be denoted at five successive points by F” 1 , F” 1 , F” 1 , 
F* ”, and F ” * . Da the starting phase, these pointB are 0, x, 2x, 3x, 
and 4x, and F Q , F^, and F” are also known. Then the five sets of 
formulas required in the starting phase are 



i s 1, 2, 3, 4 (B3) 


i » 1, 2, 3, 4 (B4) 


F^ = Fq + i XF” 



i = 1, 2, 3, 4 (B5) 
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= H 0 + ixH 




( 2 ) E .. 


i = 1, 2, 3, 4 (B6) 


x 3 

F. = F n + ixF' + ^ ' ■ x 

10 0 ^ 

D i J=0 




(B7) 


i B 1, 2, 3, 4 


(n) 


where the superscripts on the A^ and D^ 11 ^ refer to the order of 
integration. 

The constants A. and D^ 11 ) may "be read from, the following 

tables : ^ 


For 




( 1 ) 


and D 


( 1 ). 


For 




( 2 ) 


A (!) 

iJ 

»i (1) 

\J 

i\ 

0 

B 

2 

3 

4 

1 

251 

646 

-264 

106 

B 

720 

2 

29 

124 

24 

4 

B 

90 

3 

27 

102 

72 

42 

-3 

80 

4 

14 

64 

24 

64 

14 

45 

and 

V 2 ): 









(2) 



\ J 
i\ 

0 

B 

2 

3 

4 


1 

367 

540 

-282 

116 

-21 

1440 

2 

53 

144 

-30 

16 

-3 

90 

3 

441 


162 

180 

-27 

480 

4 

56 

192 

48 

64 

0 

45 
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For and : 


to 

Dl (3> 

\j 

l\ 

0 

1 

2 

3 

4 

l 

1017 

1070 

-618 

258 

-47 

10080 

2 

331 

664 

-240 

104 

-19 

630 

3 

1431 

3726 

-486 

450 

-81 

1120 

4 

744 

2176 

96 

384 

-40 

315 


It is now possible to diagram, the steps of the starting phase of 
the integration. If each bar above a function denotes an improved esti- 
mate of it, and the first estimates of FJ ' * , FI* * , F’ ' * , and F! ’ * are 

1 7 2 ’ 3 3 4 

all equal to F " * , and similarly for the H 1 1 , then the starting phase 
diagrams are u 


( 1 ). 


(F , F' , F' » F’ ' 1 , 
0 0 0 0 


F* * * , F» • * F" * F" 1 ) -* F , F' , F' 1 
1 * 2 ' 3 7 4 ' 1 * 1 * 1 


(This diagram means that the values in parentheses are used with approp- 
riate integration formulas from. (B3) to (B7) to obtain F, F', and F' ' 
at q = x. ) 


(2) . (H q , H', H”, H”, H”, H£«, H^') - 

(3) . (F 1 , Fj_, F”, B 1 , H£)h. ?£»», Ej_« 

(The preceding diagram, means that the values in parentheses are substi- 
tuted in the differential equations (Bl) and (B2) to obtain F' 1 ' and 


H 1 ' 

at q 

= X 

• ) 








(4). 


*b> 

F 6’> 

F 6’ 


F> 

", F'* 

S *i") - f 2 , 

F 2> 

F 2 

(5). 


H o j 

H 0 ' 

H- 

’ *2’ 

*b]> 

=4’ ) - 

*2’ E 2 



(6). 

< F 2’ 

F 2' 

F>', 

V 

H')-> 

F » t t 
2 





(7). 

<V 

F o> 

F 6’> 

F o' 

t F t » 

» F' 
9 2 

F ,. 

>, F”*) -> F 3 , 

F 3> 

F 3 
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DO 

O) 

-J 

O 


(8). 

(=0- 

Xb> 

TP 1 
n 0 ’ 

s r 

< S 2'- =3'. 

H4')- 

*3> *3 

(9). 


F* . 
3' 

F 1 ' 

*3 ’ 


- >f 3" 

>*h' 


do) 

< F 0> 

*b> 

F o> 

F o' 

• F' * 1 . F' 
’ 1 > 2 

"> *3 

', *i'') ^F 4 , FI 

(11) 

< E 0> 

E b> 

E o> 


’ E 3 ' * 


H 4 

(12) 


V i> 

F' 1 
■“4 > 

E 4 j 


H' ' 

* H 



It may be noted here that all four values of F' 1 1 and H' 1 have 
"been Improved, and further improvement -will require iteration of steps 
.1 to 12. The start of the second iteratibn is diagrammed as follows; 


(13) 

(14) 

(15) 

(16) 


= 0 ’ 

( ? i> 

< F 0’ F o> 


F o’ 


tt ' r 
*0 > 


F 1 1 


F" ' F* ' * , F* ' ' F' " . F* • ' ) F , F' , 

0 ’ 1 ' 2 * 3 * 4 1* 1* 

*L'> *h'> =i') ** f l> S i 

E v E[,)-* F'>t, 1“ 

F'**, F‘««, ¥$", F“>) -+ F 2 , F>, 


F 1 r 
1 


Successive sets of 12 steps are performed until the values of F^ 1 1 
and Hj,' no longer change. 


On the IBM Card Programmed Electronic Calculator, a deck of punched 
cards 2 inches thick sufficed to perform Bteps 1 to 12. Three runs of 
this starter deck at 3 minutes per run accomplished complete convergence 
in most cases . At the end of the starting process there have "been com- 


puted and stored 


F' 1 ' 
*2 > 


F“', F'", 


Hi', 


*4' 




H. , and 


=2*' 


B S' 


4 ; 
and 


S 4> 


and final estimates of F£’ * , 


H 4' 


The extension phase has now been reached. It used a different set 
of integration formulas based on the same general ideas as equations (B3) 
to (B7) . If F"', F£ M , F” ' , Fg- ' 1 , and F ” ' now designate F' 1 ' at 

any five successive points, and the subscript 5 denotes the next point, 
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F 5' = F i ,+ i tt E B j (1)j j" 

0 J=0 


H* = H! + 
5 4 


3TrE'' tl),n 


'3 3 


3=0 


*5 “ F i * +-^ T X) B J (2)jr 


3=o 

4 


t I f 

3 


H c = H, + XHJ + -fer- ^ B^ 2) HI ' 


x 5 = *4 + ^4 + "TIT 

(T ' 


5 3 


3=0 


2 3 / 3 '\ 

F 5 = F4 + XF4 + ^ F .< +-^y2^ B 3 ; 5 Mf 

C 3=0 

■where the and are given in the following table: 


BJ W 

cW 

\3 

n \ 

0 

1 

2 

3 

4 

1 

251 

-1274 

2616 

-2774 

1901 

720 

2 

135 

-692 

1446 

-1596 

1427 

1440 

3 

410 

-2116 

4476 

-5084 

5674 

20160 


The extension phase may then "be diagrammed simply as follows : 

(1) (T 4 , F£, F”, ¥' Q ", F*_”, F'”, F£««, F”*)- F g , F£, F“ 

(2) (H 4 , H«, H”, H”, H-«, H' 1 , H”)- = 5 , 

(3) (F 5 , F£, F”, H 5 , H*)- F'”, H” 


(B8) 


(B9) 


(BIO) 


(Bll) 


(B12) 
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The values of the functions at the. next point are computed, in simi- 


lar m ann er, where the latest sets of five values of I” 
used. This process advances step hy step toward infinity. 


and H' ’ are 


The extender deck of punched cards was about 3 inches thick, and 
took a little over 3 minutes per run. For Pr = 0.72, a step size of 
0.1 was used, the starting phase took 10 minutes, and the extension phase 
about 30 minutes. Men it is realized that about '11,000 operations were 
performed in the 40 minutes per run, it may be seen that solution of the 
present problem would have been prohibitively difficult on desk-type cal- 
culators . Simplifications in method would have sacrificed accuracy or 
required smaller step size. 

In two-point boundary value problems where one point is infinity, 
some problems of judgment are involved as to where infinity is, and as to 
when a satisfactory approximation to a solution has been obtained. In 
most cases this question was settled for the present problem by calling 
a run satisfactory when it fell between two runs for which F' and H 
did not differ at important points in the fourth decimal place, and for 
which F' and H flattened out at zero, correct to four decimal places. 

Certain difficulties were met in the attempt to use the basic pro- 
cedure previously discussed. These necessitated certain modifications. 

For Pr = 2, 10, 100, and 1000, H would settle down to zero at an 
early stage, but while F 1 was still coming down, H 1 * would begin to 
oscillate and these oscillations increased and fed back into all other 
functions.. This trouble was avoided by the following modifications: It 

is a consequence of equation (B2) that 


H’(q) = H'(0)exp^3 Pr F(t) dtj 


H* (tj - x )exp ^-3 


Pr F(t) dtj 


(B13) 


The extension phase was modified to require the additional integra- 
tion formulas 

5 



F(t) dt = — 


720 




(B14) 


i=l 



J 

. 1 
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I 


5 

% = =4 + ffo Y, A i (B15) 

1=1 

where = -19, Ag = 106, A^ = -264, A^_ = 646, and A^ = 251. 

These fCrmulas were used along with equations (B8) to (BIO) according 
to the following diagram; 


(1) 

0V 



JM 1 
0 

t 

y 

F" * 

, F" 

(2) 

(F 1' 


F 3> 

F 4> 

F 5 ; 

- Hp 

- H 5 

( 3 ) 

(H', 


* 3 ’ 

H 4> 


. V 

** H 5 

( 4 ) 

<*5' 

F 5' 

* 5 ’ 

H 5 ) 

-* 

JM 1 t 

5 



The 

value F^ 

t t 

is 

discarded 


, F”*, F"«) - F 5 , Fj., F“ 

"by means of (B14) and (B13) 

"by means of (B15) 

a-nri F 1 1 1 at the last five points is 


used to repeat the whole process again and again ad infinitum. Ab long 
as F stays positive, H' is guaranteed to approach zero and H will 
flatten out to some value and not oscillate. 


For Er = 0.01, 0.72, 22/30, and 1, the F* 1 ’ Began to oscillate 
at an advanced point and these oscillations grew and fed into the other 
functions. For all cases hut Pr = 0.01, the oscillations appeared very 
late, near the end of the run, and a suitable halving of step size when 
oscillation was detected in the fourth differences of F ,,f was suffic- 
ient to avoid the difficulty. But for the 0.01 case, oscillations of 
F 1 1 1 appeared early in the run, namely, soon after the peak in F. 

These oscillations were found to he step-size connected, so that reduc- 
tion of the step to 0.02 avoided them. Even then oscillations in F MI 
would begin to appear e^ery 25 steps or so, and these were smoothed out 
regularly by repeated r uns of a deck similar to the starter deck. Each 
run under these conditions took about 16 hours, making this the most 
difficult case to solve. 


< 
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APPENDIX 0 

DERIVATION OF FLOW PARAMETER 


By definition the shear stress is given By 

T = "°© 0 . 


(Cl) 


To express (Sij/Sy)q in terms of the known function F(t]), use can he 
made of equations (60) and (62). Then 


■Su 

dY 




Substitution of this expression into equation (Cl) yields the flow 
parameter 


=F» '(0) 

( 4 Gr x 3 ) 4 ^- 'V* 2 ) 

Note that from, the general derivation, the flow parameter contains the 
viscosity evaluated at two different points. Becall, however, that the 
analysis has shown that to a first approximation the variation of vis- 
cosity with temperature can he neglected. Thus the Viscosity can he 
taken as constant in the entire flow field. 
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APPENDIX D 

DERIVATION OF BEAT -TRANSFER PARAMETER 


The local Nusselt number is defined as 


s hi _ -X /St\ 
k ” (T 0 -To) \ST/ o 


(Dl) 


To express (St/Sy)^ 
of equations (6l) and 


in terms of the known function 


(62). 

Thus 

Bt _ 

(To-t.) 

Sy “ 

X 



H(tj) , use is made 


Substitution of this expression into equation (Dl) yields the heat- 
transfer parameter 


— y = - H> (0) (D2) 

(Gr z /4) 4 

The heat -transfer parameter as given "by equation (D2) is, as was pre- 
viously stated, a local parameter. It is often desired to compute the 
average (over the length X) value of this parameter. To this end, the 
Nusselt number (as given in equation (Dl)) must he defined in terms of 
an average heat -transfer coefficient and the quantity thus obtained must 
then he integrated over the length X and divided hy X. This procedure 
yields the result 

nu . | (nu) aT 

It is from this last equation that the 0.75 reduction factor previously 
discussed was obtained. 
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IABIE I - FEffiCTIOIS 


F AND H AND DERIVATIVES FOB VARIOUS ERAKDTL NUMBERS 


(a) Brandt! number, O.Ol 


n 

F 

F ' 

F ' ' 

H 

H ' 


F 

F ' 

JP * * 

H 

H ' 

0 

0.0000 

0.0000 

0.9862 

1.0000 

- 0.0812 

4.16 

1.9319 

0.4037 

- 0.0615 

0.6741 

- 0.0723 

.1 

.0048 

.0936 

.8868 

.9919 

-.0812 

4.18 

1.9400 

.4024 

-.0614 

.6726 

-.0722 

.2 

.0184 

.1774 

.7891 

.9838 

-.0812 

4.20 

1.9480 

.4012 

-.0612 

.6712 

-.0721 

.3 

.0399 

.2516 

.6942 

.9756 

-.0812 

4.24 

1.9640 

.3988 

-.0609 

.6683 

-.0720 

.4 

.0684 

.3164 

.6030 

.9675 

-.0811 

4.28 

1.9799 

.3964 

-.0605 

.6654 

-.0718 

.5 

.1029 

.3723 

.5163 

.9594 

-.0811 

4.32 

1.9957 

.3939 

-.0602 

.6626 

-.0716 

.6 

.1426 

.4198 

.4349 

.9513 

-.0811 

4.36 

2.0114 

.3915 

-.0599 

.6597 

-.0714 

.7 

.1866 

.4595 

.3595 

.9432 

-.0811 

4.40 

2.0270 

.3891 

-.0595 

.6568 

-.0713 

.8 

.2342 

.4919 

.2906 

.9351 

-.0810 

4.50 

2.0657 

.3832 

-.0587 

.6497 

-.0708 

.9 

.2848 

.5178 

.2285 

.9270 

-.0809 

4.60 

2.1037 

.3774 

-.0579 

.6427 

-.0704 

1.0 

.3376 

.5379 

.1734 

.9189 

-.0809 

4.70 

2.1411 

.3716 

-.0572 

.6357 

-.0699 

1.1 

.3922 

.5527 

.1253 

.9108 

-.0808 

4.80 

2.1780 

.3660 

-.0564 

.6287 

-.0695 

1.2 

.4480 

.5631 

.0839 

.9028 

-.0807 

4.90 

2.2143 

.3604 

-.0556 

.6218 

-.0690 

1.3 

.5047 

.5697 

.0489 

.8947 

-.0806 

5.00 

2.2501 

.3548 

-.0549 

.6149 

-.0686 

1.4 

.5618 

.5731 

.0198 

.8866 

-.0804 

5.10 

2.2853 

.3494 

-.0541 

.6081 

-.0681 

1.5 

.6192 

.5739 

- 1.0039 

.8786 

-.0803 

5.20 

2.3200 

.3440 

-.0534 

.6013 

-.0676 

1.6 

.6765 

.5725 

-.0229 

.8706 

-.0801 

5.40 

2.3877 

.3335 

-.0520 

.5878 

-.0667 

1.7 

.7336 

.5694 

-.0379 

.8626 

-.0800 

5.60 

2.4534 

.3232 

-.0506 

.5746 

-.0657 

1.8 

.7904 

.5650 

-.0493 

.8546 

-.0798 

5.80 

2.5170 

.3132 

-.0493 

.5615 

-.0648 

1.9 

.8466 

.5596 

-.0579 

.8466 

-.0796 

6.00 

2.5787 

.3035 

-.0479 

.5487 

-.0638 

2.0 

.9023 

.5535 

-.0641 

.8387 

-.0794 

6.20 

2.6384 

.2940 

-.0467 

.5360 

-.0628 

2.1 

.9573 

.5469 

-.0685 

.8307 

-.0792 

6.40 

2.6963 

.2849 

-.0454 

.5236 

-.0618 

2.2 

1.0117 

.5399 

-.0714 

.8228 

-.0789 

6.60 

2.7524 

.2759 

-.0442 

.5113 

-.0608 

2.3 

1.0653 

.5326 

-.0732 

.8150 

-.0787 

6.80 

2.8067 

.2671 

-.0430 

.4993 

-.0598 

2.4 

1.1182 

.5253 

-.0742 

.8071 

-.0784 

7.00 

2.8593 

.2586 

-.0419 

.4874 

-.0588 

2.5 

1.1703 

.5178 

-.0746 

.7993 

-.0782 

7.40 

2.9594 

.2423 

-.0397 

.4643 

-.0568 

2.6 

1.2217 

.5104 

-.0745 

.7915 

-.0779 

7.80 

3.0532 

.2269 

-.0376 

.4420 

-.0547 

2.7 

1.2724 

.5029 

-.0741 

.7837 

-.0776 

8.20 

3.1411 

.2123 

-.0356 

.4205 

-.0527 

2.8 

1.3223 

.4956 

-.0735 

.7760 

-.0773 

8.60 

3.2232 

.1984 

-.0337 

.3998 

-.0508 

2.9 

1.3715 

.4882 

-.0728 

.7682 

-.0770 

9.00 

3.2999 

.1853 

-.0319 

.3799 

-.0489 

3.0 

1.4200 

.4810 

-.0719 

.7606 

-.0766 

9.40 

3.3715 

.1729 

-.0302 

.3607 

-.0470 

3.1 

1.4677 

.4739 

-.0711 

.7529 

-.0763 

9.80 

3.4383 

.1612 

-.0285 

.3423 

-.0451 

3.2 

1.5148 

.4668 

-.0701 

.7453 

-.0760 

10.20 

3.5005 

.1501 

-.0270 

.3246 

-.0432 

3.3 

1.5611 

.4598 

-.0692 

.7377 

-.0756 

10.60 

3.5584 

.1396 

-.0254 

.3069 

-.0413 

3.4 

1.6067 

.4530 

-.0683 

.7302 

-.0753 

11.00 

3.6123 

.1297 

-.0241 

.2915 

-.0397 

3.45 

1.6293 

.4496 

-.0678 

.7264 

-.0751 

11.40 

3.6622 

.1203 

-.0228 

.2759 

-.0380 

3.5 

1.6517 

.4462 

-.0674 

.7227 

-.0749 

11.80 

3.7086 

.1114 

-.0215 

.2610 

-.0364 

3.55 

1.6739 

.4428 

-.0669 

.7189 

-.01 in 

12.20 

3.7514 

.1031 

-.0203 

.2468 

-.0348 

3.6 

1.6960 

.4395 

-.0664 

.7152 

-.0745 

12.60 

3.7911 

.0952 

-.0192 

.2332 

-.0332 

3.65 

1.7179 

.4362 

-.0660 

.7115 

-.0743 

13.00 

3.8276 

.0877 

->.0181 

.2202 

-.0317 

3.7 

1.7396 

.4329 

-.0655 

.7078 

-.0741 

13.60 

3.8771 

.0773 

-.0165 

.2018 

-.0296 

3.75 

1.7611 

.4296 

-.0651 

.7041 

-.0739 

14.20 

3.9206 

.0679 

-.0151 

.1847 

-.0276 

3.8 

1.7825 

.4264 

-.0646 

.7004 

-.0738 

14.80 

3.9587 

.0592 

-.0138 

.1687 

-.0257 

3.85 

1.8308 

.4232 

-.0642 

.6967 

-.0736 

15.40 

3.9918 

.0513 

-.0126 

.1538 

-.0239 

3.9 

1.8249 

.4200 

-.0638 

.6930 

-.0734 

16.00 

4.0204 

.0441 

-.0114 

.1399 

-.0223 

3.95 

1.8458 

.4168 

-.0633 

.6893 

-.0732 

17.00 

4.0591 

.0335 

-.0097 

.1189 

-.0197 

4.00 

1.8665 

.4136 

-.0629 

.6857 

-.0729 

18.00 

4.0880 

.0246 

-.0082 

.1003 

-.0175 

4.05 

1.8871 

.4105 

-.0625 

.6821 

-.0727 

19.00 

4.1088 

.0171 

-.0069 

.0839 

-.0154 

4.10 

1.9076 

.4074 

-.0620 

.6784 

-.0725 

20.00 

4.1226 

.0108 

-.0057 

.0694 

-.0137 

4.12 

1.9157 

.4061 

-.0619 

.6770 

-.0725 

21.00 

4.1308 

.0057 

-.0046 

.0565 

-.0121 

4.14 

1.9238 

.4049 

-.0617 

.6755 

-.0724 

22.00 

4.1343 

.0015 

-.0037 

.0452 

-.0107 
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TABLE I - FOHCTIOHS F AUD H ADD DERIVATIVES 


(t) FrsncLtl number, 0.72. (oj Prandtl number, 0.733 


n 

F 

F' 

F » * 

H 

H ’ 


_J1_ 

F 

F' 

F' ' 

E 

H' 

0 

0.0000 

0.0000 

0.6760 

1.0000 

- 0.5046 


0 

0.0000 

0.0000 

0.6741 

1.0000 

- 0.5080 

.1 

.0032 

.0627 

.5785 

.9495 

-.5045 


.1 

.0032 

.0625 

.5767 

.9492 

-.5079 

.2 

.0122 

.1159 

.4866 

.8991 

-.5037 


.2 

.0122 

.1155 

.4849 

.8984 

-.5071 

.3 

.0261 

.1602 

.4007 

.8488 

-.5016 


.3 

.0260 

.1597 

.3990 

.8478 

-.5050 

.4 

.0440 

.1962 

.3210 

.7989 

-.4979 


.4 

.0438 

.1955 

.3194 

.7975 

-.5012 

.5 

.0651 

.2246 

.2479 

.7493 

-.4921 


.5 

^0649 

.2238 

.2465 

.7477 

-.4953 

.6 

.0887 

.2460 

.1817 

.7005 

-.4840 


.6 

.0884 

.2451 

.1804 

.6985 

-.4870 

.7 

.1141 

.2612 

.1224 

.6526 

-.4735 


.7 

.1137 

.2601 

.1212 

.6503 

-.4763 

.8 

.1407 

.2708 

.0700 

.6059 

-.4607 


.8 

.1402 

.2695 

.0691 

.6033 

-.4632 

.9 

.1681 

.2754 

.0246 

.5606 

-.4456 


.9 

.1674 

.2741 

.0237 

.5578 

-.4478 

1.0 

.1957 

.2759 

-.0143 

.5168 

-.4284 


1.0 

.1949 

.2745 

-.0149 

.5139 

‘ -.4303 

1.1 

.2231 

.2728 

-.0468 

.4749 

-.4095 


1.1 

.2222 

.2713 

-.0473 

.4718 

-.4110 

1.2 

.2501 

.2667 

-.0734 

.4350 

-.3891 


1.2 

.2490 

.2652 

-.0737 

.4317 

-.3902 

1.3 

.2764 

.2583 

-.0945 

.3972 

-.3676 


1.3 

.2752 

.2568 

-.0946 

.3938 

-.3684 

1.4 

.3017 

.2480 

-.1106 

.3615 

-.3453 


1.4 

.3003 

.2465 

-.1106 

.3581 

-.3458 

1.5 

.3260 

.2363 

-.1224 

.3281 

-.3227 


1.5 

.3244 

.2348 

-.1222 

.3246 

-.3228 

1.6 

.3490 

.2236 

-.1302 

.2970 

-.3000 


1.6 

.3473 

.2222 

-.1299 

.2935 

-.2998 

1.7 

.3707 

.2104 

-.1347 

.2681 

-.2775 


1.7 

.3688 

.2090 

-.1342 

.2647 

-.2771 

1.8 

.3910 

.1968 

-.1363 

.2415 

-.2556 


1.8 

.3891 

.1954 

-.1358 

.2381 

-.2549 

1.9 

.4100 

.1832 

-.1356 

.2170 

-.2344 


1.9 

.4079 

'.1819 

-.1350 

.2136 

-.2335 

2.0 

.4277 

.1697 

-.1331 

.1945 

-.2141 


2.0 

.4254 

.1685 

-.1324 

.1913 

-.2130 

2.1 

.4440 

.1566 

-.1290 

.1741 

-.1949 


2.1 

.4416 

.1554 

-.1283 

.1710 

-.1937 

2.2 

.4590 

.1440 

-.1238 

.1555 

-.1768 


2.2 

.4565 

.1429 

-.1230 

.1526 

-.1754 

2.3 

.4728 

.1319 

-.1178 

.1387 

-.1598 


2.3 

.4702 

.1309 

-.1170 

.1359 

-.1584 

2.4 

.4854 

.1204 

-.1113 

.1235 

-.1441 


2.4 

.4827 

.1195 

-.1104 

.1208 

-.1427 

2.5 

.4969 

.1097 

-.1044 

.1098 

-.1296 


2.5 

.4941 

.1088 

-.1035 

.1073 

-.1281 

. 2.6 

.5074 

.0996 

-.0973 

.0976 

-.1163 


2.6 

.5045 

.0988 

-.0965 

.0952 

-.1148 

2.7 

.5168 

.0902 

-.0903 

.0865 

-.1041 


2.7 

.5139 

.0895 

-.0895 

.0843 

-.1026 

2.8 

.5254 

.0815 

-.0834 

.0767 

-.0930 


2.8 

.5224 

.0809 

-.0826 

.0746 

-.0916 

2.9 

.5332 

.0735 

-.0767 

.0679 

-.0830 


2.9 

.5301 

.0729 

-.0759 

.0660 

-.0816 

3.0 

.5401 

.0661 

-.0703 

.0601 

-.0739 


3.0 

.5370 

.0657 

-.0695 

.0583 

-.0725 

3.1 

.5464 

.0594 

-.0642 

.0531 

-.0657 


3.1 

.5433 

.0590 

-.0635 

.0514 

-.0644 

3.2 

.5520 

.0533 

-.0585 

.0469 

-.0584 


3.2 

.5489 

.0530 

-.0578 

.0454 

-.0571 

3.3 

.5571 

.0477 

-.0531 

.0414 

-.0518 


3.3 

.5539 

.0475 

-.0524 

.0400 

-.0506 

3.4 

.5616 

.0427 

-.0481 

.0365 

-.0459 


3.4 

.5584 

.0425 

-.0474 

.0352 

-.0448 

3.6 

.5692 

.0339 

-.0392 

.0284 

-.0359 


3.6 

.5660 

.0339 

-.0386 

.0273 

-.0350 

3.8 

.5735 

.0269 

-.0317 

.0220 

-.0281 


3.8 

.5720 

.0269 

-.0312 

.0211 

-.0272 

4.0 

.5801 

.0212 

-.0255 

.0170 

-.0219 


4.0 

.5769 

.0213 

-.0250 

.0163 

-.0211 

4.2 

.5838 

.0166 

-.0204 

.0132 

-.0170 


4.2 

.5807 

.0169 

-.0200 

.0126 

-.0164 

4.4 

.5868 

.0130 

-.0162 

..0102 

-.0132 


4.4 

.5837 

.0133 

-.0159 

.0097 

-.0127 

4.6 

.5891 

.0101 

-.0128 

.0078 

-.0102 


4.6 

.5860 

.0105 

-.0125 

.0074 

-.0098 

4.8 

.5908 

.0078 

-.0102 

.0060 

-.0080 


4.8 

.5879 

.0082 

-.0099 

.0057 

-.0076 

5.0 

.5922 

.0060 

-.0079 

.0046 

-.0061 


5.0 

.5893 

.0065 

-.0077 

.0044 

-.0058 

5.2 

.5933 

.0046 

-.0063 

.0035 

-.0048 


5.2 

.5905 

.0051 

-.0061 

.0033 

-.0045 

5.4 

.5941 

.0034 

-.0048 

.0027 

-.0037 


5.4 

.5914 

.0040 

-.0047 

.0025 

-.0035 

5.8 

.5951 

.0019 

-.0029 

.0015 

-.0022 


5.6 

.5921 

.0032 

-.0036 

.0019 

-.0027 

6.2 

.5957 

.0010 

-.0018 

.0009 

-.0013 


5.8 

.5927 

.0025 

-.0028 

.0015 

-.0021 

6.8 

.5960 

.0003 

-.0008 

.0003 

-.0006 


6.0 

.5932 

.0021 

-.0022 

.0011 

-.0016 

7.3 

.5961 

.0000 

-.0004 

.0001 

-.0003 


6.4 

.5938 

.0014 

-.0013 

.0006 

-.0009 








6.8 

.5943 

.0010 

-.0008 

.0003 

-.0006 








7.2 

.5946 

.0007 

-.0004 

.0001 

-.0003 








7.6 

.5949 

.0006 

-.0002 

. 0000 ' 

-.0002 








8.0 

.5951 

.0005 

-.0001 

.0000 

-.0001 
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FOE VABI0U3 ERAHDTL HUMBSRS - CQtfi .' JJilM ) 


( 4 ) Prandtl number , 1 


n 

F 

F » 

F ' 1 

H 

H ' 

0 

0.0000 

0.0000 

0.6421 

1.0000 

- 0.5671 

.1 

.0030 

.0593 

.5450 

.9433 

-.5669 

.2 

.0115 

.1092 

.4540 

.8867 

-.5657 

.3 

.0246 

.1503 

.3694 

.8302 

-.5627 

.4 

.0413 

.1833 

.2916 

.7742 

-.5572 

.5 

.0610 

.2089 

.2208 

.7189 

-.5488 

.6 

.0829 

.2277 

.1572 

.6645 

-.5371 

.7 

.1064 

.2406 

.1008 

.6116 

-.5221 

.8 

.1308 

.2481 

.0516 

.5602 

-.5038 

.9 

.1558 

.2511 

.0093 

.5109 

-.4826 

1.0 

.1809 

.2502 

-.0263 

.4638 

-.4589 

1.1 

.2058 

.2461 

-.0557 

.4192 

-.4330 

1.2 

.2300 

.2393 

-.0793 

.3772 

-.4056 

1.3 

.2535 

.2304 

-.0975 

.3381 

-.3772 

1.4 

.2761 

.2199 

-.1110 

.3018 

-.3484 

1.5 

.2975 

.2083 

-.1203 

.2684 

-.3197 

1.6 

.3177 

.1960 

-.1260 

.2379 

-.2915 

1.7 

.3367 

.1832 

-.1287 

.2101 

-.2642 

1.8 

.3543 

.1703 

-.1288 

.1850 

-.2382 

1.9 

.3707 

.1575 

-.1268 

.1624 

-.2136 

2.0 

.3859 

.1450 

-.1233 

.1422 

-.1907 

2.1 

.3997 

.1329 

-.1185 

.1242 

-.1695 

2.2 

.4125 

.1213 

-.1127 

.1082 

-.1501 

2.3 

.4240 

.1104 

-.1064 

.0941 

-.1324 

2.4 

.4346 

.1001 

-.0997 

.0817 

-.1164 

2.5 

.4441 

.0904 

-.0928 

.0708 

-.1020 

2.6 

.4527 

.0815 

-.0859 

.0613 

-.0892 

2.7 

.4604 

.0733 

-.0791 

.0529 

-.0777 

2.8 

.4673 

.0657 

-.0725 

.0457 

-.0676 

2.9 

.4736 

.0588 

-.0662 

.0392 

-.0587 

3.0 

.4791 

.0524 

-.0602 

.0339 

-.0509 

3.1 

.4841 

.0467 

-.0546 

.0291 

-.0441 

3.2 

.4885 

.0415 

-.0493 

.0250 

-.0381 

3.3 

.4924 

.0368 

-.0444 

.0215 

-.0329 

3.4 

.4959 

.0326 

-.0399 

.0185 

-.0283 

3.6 

.5016 

.0254 

-.0321 

.0136 

-.0210 

3.8 

.5061 

.0197 

-.0255 

.0099 

-.0155 

4.0 

.5096 

.0151 

-.0202 

.0072 

-.0115 

4.2 

.5122 

.0116 

-.0158 

.0053 

-.0084 

4.4 

.5143 

.0087 

-.0124 

.0038 

-.0062 

4.6 

.5158 

.0066 

-.0096 

.0027 

-.0045 

4.8 

.5169 

.0049 

-.0075 

.0020 

-.0034 

5.0 

.5177 

.0035 

-.0057 

.0014 

-.0024 

5.2 

.5183 

.0025 

-.0044 

.0010 

-.0018 

5.5 

.5189 

.0014 

-.0029 

.0006 

-.0011 

6.0 

.5194 

.0004 

-.0014 

.0002 

-.0005 

6.25 

.5194 

.0000 

-.0010 

.0000 

-.0004 


( e ) Prandtl number , 2 


in 

F 

F ' 

F ' • 

H 


0 

0.0000 

0 .0000 

0.5713 

1.0000 

- 0.7165 

.1 

.0027 

.0523 

.4749 

.9284 

-.7161 

.2 

.0101 

.0952 

.3861 

.8569 

-.7135 

.3 

.0215 

.1297 

.3049 

.7858 

-.7069 

.4 

.0358 

.1565 

.2318 

.7157 

-.6949 

.5 

.0525 

.1763 

.1866 

.6470 

-.6768 

.6 

.0709 

.1901 

.1095 

.5805 

-.6523 

.7 

.0904 

.1985 

.0602 

.5168 

-.6215 

.8 

.1104 

.2024 

.0185 

.4564 

-.5852 

.9 

.1307 

.2024 

-.0161 

.3999 

-.5443 

1.0 

.1508 

.1994 

-.0440 

.3476 

-.5002 

1.1 

.1705 

.1938 

-.0659 

.2999 

-.4543 

1.2 

.1895 

.1864 

-.0824 

.2568 

-.4077 

1.3 

.2077 

.1775 

-.0942 

.2183 

-.3619 

1.4 

.2250 

.1677 

-.1020 

.1844 

-.3178 

1.5 

.2412 

.1572 

-.1063 

.1547 

-.2763 

1.6 

.2564 

.1465 

-.1078 

.1290 

-.2380 

1.7 

.2706 

.1357 

-.1071 

.1069 

-.2032 

1.8 

.2836 

.1252 

-.1046 

.0882 

-.1721 

1.9 

.2956 

.1149 

-.1008 

.0724 

-.1446 

2.0 

.3066 

.1050 

-.0960 

.0592 

-.1207 

2.1 

.3166 

.0957 

-.0906 

.0482 

-.1001 

2.2 

.3257 

.0869 

-.0849 

.0391 

-.0826 

2.3 

.3340 

.0787 

-.0789 

.0316 

-.0677 

2.4 

.3415 

.0711 

-.0729 

.0254 

-.0553 

2.5 

.3483 

.0641 

-.0671 

.0204 

-.0450 

2.6 

.3543 

.0577 

-.0614 

.0164 

-.0364 

2.7 

.3598 

.0518 

-.0560 

.0131 

-.0294 

2.8 

.3647 

.0465 

-.0509 

.0105 

-.0237 

2.9 

.3691 

.0417 

-.0461 

.0083 

-.0190 

3.0 

.3731 

.0373 

-.0416 

.0066 

-.0152 

3.1 

.3766 

.0333 

-.0375 

.0053 

-.0121 

3.2 

.3798 

.0298 

-.0338 

.0042 

-.0097 

3.3 

.3826 

.0266 

-.0303 

.0033 

-.0077 

3.4 

.3851 

.0237 

-.0272 

.0026 

-.0061 

3.6 

.3893 

.0188 

-.0218 

.0017 

-.0038 

3.8 

.3927 

.0149 

-.0174 

.0010 

-.0024 

4.0 

.3953 

.0118 

-.0138 

.0007 

-.0015 

4.2 

.3974 

.0094 

-.0109 

.0004 

-.0009 

4.4 

.3991 

.0074 

-.0086 

.0003 

-.0006 

4.6 

.4004 

.0059 

-.0068 

.0002 

-.0004 

4.8 

.4015 

.0047 

-.0054 

.0001 

-.0002 

5.0 

.4023 

.0037 

-.0042 

.0001 

-.0001 

5.4 

.4035 

.0024 

-.0026 

.0001 

-.0001 

5.8 

.4043 

.0015 

-.0016 

.0000 

.0000 

6.4 

.4049 

.0008 

-.0008 

.0000 

.0000 

7.0 

.4053 

.0005 

-.0004 

.0000 

.0000 

8.0 

.4056 

.0002 

-.0001 

.0000 

.0000 

9.0 

.4058 

.0001 

,0000 

.0000 

.0000 

10.0 

.4059 

.0001 

.0000 

.0000 

.0000 

11.0 

.4059 

.0001 

.0000 

.0000 

.0000 
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TABLE I - FUHCTIOHS F 


ARD H AJtD DERI7ATTVB3 


(f 

Prandtl number , 10 

A - 

F 

F ' 

Jit t 

H 

H ' 

0 

0 .0000 

0.0000 

0.4192 

1.0000 

- 1.1694 

.1 

.0019 

.0371 

.3251 

.8831 

- 1.1671 

.2 

.0071 

.0654 

.2428 

.7670 

- 1.1521 

.3 

.0147 

.0861 

.1723 

.6534 

- 1.1155 

.4 

.0241 

.1003 

.1134 

.5448 

- 1.0526 

.5 

.0346 

.1091 

.0655 

.4437 

-.9640 

.6 

.0458 

.1137 

.0279 

.3527 

-.8545 

.7 

.0573 

.1150 

-.0011 

.2733 

-.7322 

.8 

.0687 

.1137 

-.0221 

.2064 

-.6061 

.9 

.0800 

.1107 

-.0367 

.1519 

-.4849 

1.0 

.0908 

.1066 

-.0462 

.1090 

-.3753 

1.1 

.1012 

.1016 

-.0518 

.0763 

-.2813 

1.2 

.1111 

.0963 

-.0545 

.0522 

-.2045 

1.3 

.1205 

.0908 

-.0552 

.0349 

-.1445 

1.4 

.1293 

.0853 

-.0544 

.0228 

-.0993 

1.5 

.1376 

.0799 

-.0527 

.0146 

-.0665 

1.6 

.1453 

.0748 

-.0505 

.0092 

-.0435 

1.7 

.1525 

.0699 

-.0480 

.0056 

-.0278 

1.8 

.1593 

.0652 

-.0453 

.0034 

-.0174 

1.9 

.1656 

.0608 

-.0427 

.0020 

-.0107 

2.0 

.1714 

.0567 

-.0400 

.0012 

-.0065 

2.1 

.1769 

.0528 

-.0375 

.0007 

-.0038 

2.2 

.1820 

.0491 

-.0351 

.0004 

-.0022 

2.3 

.1868 

.0458 

-.0328 

.0002 

-.0013 

2.4 

.1912 

.0426 

-.0306 

.0001 

-.0007 

2.5 

.1953 

.0396 

-.0286 

.0001 

-.0004 

2.6 

.1991 

.0369 

-.0267 

.0000 

-.0002 

2.7 

.2027 

.0343 

-.0249 

.0000 

-.0001 

2.8 

.2060 

.0319 

-.0232 

.0000 

-.0001 

2.9 

.2090 

.0297 

-.0216 

.0000 

.0000 

3.0 

.2119 

.0276 

-.0201 

.0000 

.0000 

3.1 

.2146 

.0256 

-.0187 

.0000 

.0000 

3.2 

.2170 

.0238 

-.0174 

.0000 

.0000 

3.3 

.2193 

.0221 

-.0162 

.0000 

.0000 

3.4 

.2215 

.0206 

-.0151 

.0000 

.0000 

3.6 

.2253 

.0178 

-.0131 

.0000 

.0000 

3.8 

.2286 

.0153 

-.0113 

.0000 

.0000 

4.0 

.2314 

.0132 

-.0098 

.0000 

.0000 

4.2 

.2339 

.0114 

-.0084 

.0000 

.0000 

4.4 

.2360 

.0098 

-.0073 

.0000 

.0000 

4.6 

.2379 

.0085 

-.0063 

.0000 

.0000 

4.8 

.2394 

.0073 

-.0054 

.0000 

.0000 

5.0 

• 2 - J 08 

.0063 

-.0047 

.0000 

.0000 

5.4 

.2430 

.0047 

-.0035 

.0000 

. 0000 ' 

5.8 

.2446 

.0035 

-.0026 

.0000 

.0000 

6.2 

.2458 

.0026 

-.0019 

.0000 

.0000 

7.0 

.2474 

.0014 

-.0011 

.0000 

.0000 

8.0 

.2484 

.0007 

-.0005 

.0000 

.0000 

9.0 

.2489 

.0003 

-.0002 

..0000 

.0000 

10.0 

.2491 

.0002 

-.0001 

.0000 

.0000 


( g ) Prandtl number , 100 


1 ) 

F 

F ' 

F 1 1 

H 

H 1 

0.000 

0 .000 

0 .0000 

0.2517 

1.0000 

- 2.191 

.025 

.0001 

.0060 

.2274 

. 9452 - 

- 2. 191 

.050 

.0003 

.0114 

.2044 

.8905 

- 2.188 

.075 

.0006 

.0162 

.1828 

.8359 

- 2.180 

.100 

.0011 

.0205 

.1626 

.7815 

- 2.166 

.125 

.0017 

.0244 

.1438 

.7276 

- 2.144 

.150 

.0023 

.0277 

.1263 

.6744 

- 2.113 

.175 

.0031 

.0307 

.1101 

.6221 

- 2.071 

.200 

.0039 

.0332 

.0952 

.5709 

- 2.018 

.225 

.0047 

.0355 

.0816 

.5213 

- 1.954 

.250 

.0056 

.0373 

.0692 

.4733 

- 1 ,880 

.275 

.0066 

.0389 

.0580 

.4273 

- 1.786 

.300 

.0076 

.0402 

.0479 

.3836 

- 1.704 

.325 

.0086 

.0413 

.0389 

.3422 

- 1.604 

.350 

.0096 

.0422 

.0309 

.3034 

- 1.498 

.375 

.0107 

.0429 

.0238 

.2674 

- 1.388 

.400 

.0118 

.0434 

.0176 

.2341 

- 1.276 

.425 

.0129 

.0438 

.0123 

.2036 

- 1.163 

.450 

.0140 

.0440 

.0076 

.1759 

- 1.052 

.475 

.0151 

. 0442 ' 

.0036 

.1510 

-.9430 

.500 

.0162 

.0442 

.0002 

.1287 

-.8393 

.525 

.0173 

.0442 

-.0026 

.1089 

-.7404 

.550 

.0184 

.0441 

-.0050 

.0916 

-.6478 

.575 

.0195 

.0439 

-.0070 

.0765 

-.5621 

.600 

.0206 

.0437 

-.0087 

.0634 

-.4837 

.65 

.0227 

.0432 

-.0111 

.0427 

-.3496 

.70 

.0249 

.0426 

-.0126 

.0280 

-.2440 

.75 

.0270 

.0420 

-.0135 

.0178 

-.1657 

.80 

.0291 

.0413 

-.0140 

.0111 

-.1088 

.85 

.0311 

.0406 

-.0142 

.0087 

-.0693 

.90 

.0331 

.0399 

-.0142 

.0039 

-.0428 

.95 

.0351 

.0392 

-.0141 

.0022 

-.0256 

1.00 

.0371 

' .0385 

-.0140 

.0012 

-.0149 

1.10 

.0408 

.0371 

-.0136 

.0004 

-.0046 

1.20 

.0445 

.0357 

-.0132 

.0001 

-.0013 

1.30 

.0480 

.0344 

-.0128 

.0000 

-.0003 

1.40 

.0514 

.0332 

-.0124 

.0000 

-.0001 

1.50 

.0546 

.0320 

-.0119 

.0000 

.0000 

1.60 

.0578 

.0308 

-.0116 

.0000 

.0000 

1.70 

.0608 

.0297 

-.0112 

.0000 

.0000 

1.80 

.0637 

.0286 

-.0108 

.0000 

.0000 

1.90 

.0665 

.0275 

-.0104 

.0000 

.0000 

2.0 

.0692 

.0265 

-.0101 

.0000 

.0000 

2.1 

.0718 

.0255 

-.0097 

.0000 

.0000 

2.2 

.0743 

.0245 

-.0094 

.0000 

.0000 

2.3 

.0767 

.0236 

-.0091 

.0000 

.0000 

2.4 

.0790 

.0227 

-.0088 

.0000 

.0000 

2.6 

.0834 

.0210 

-.0082 

.0000 

.0000 

2.8 

.0874 

.0195 

-.0076 

.0000 

.0000 

3.0 

.0912 

.0180 

-.0071 

.0000 

.0000 
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FOB VARIOUS ERAJJDEL NUMBERS - CONCLUDED 


(g) Pran&tl number, 100 - Conolu&ed 


n 

F 

F ' 

» 

H 

H * 

3.2 

3.4 

3.6 

3.8 
4 . 0 - 

4.4 

4.8 

5.2 

5.6 
6.0 

6.6 

7.2 
8.0 
9.0 

10.0 

11.0 

12.0 

13.0 

0.0947 

.0979 

.1008 

.1035 

.1061 

.1105 

.1143 

.1176 

.1203 

.1226 

.1254 

.1276 

.1297 

.1315 

.1326 

.1332 

.1335 

.1336 

0.0166 

.0154 

.0142 

.0131 

.0121 

.0103 

.0088 

.0074 

.0063 

.0053 

.0041 

.0032 

.0022 

.0014 

.0008 

.0004 

.0002 

.0000 

- 0.0066 

-.0061 

-.0056 

-.0052 

-.0049 

-.0042 

-.0036 

-.0031 

-.0026 

-.0022 

-.0018 

-.0014 

-.0010 

-.0007 

-.0005 

-.0003 

-.0002 

-.0001 

0.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

0.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

,.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 


(i) Prartdtl number, 1000 



F 

F ’ 

F ' • 

H 

H ' 

0. 

0.0000 

0.0000 

0.1450 

1.0000 

- 3.966 

.025 

.0000 

.0033 

.1212 

.9009 

- 3.962 

.050 

.0002 

.0061 

.0999 

.8021 

- 3.933 

.075 

.0003 

.0083 

.0811 

.7046 

- 3.861 

.100 

.0006 

.0102 

.0647 

.6096 

- 3.731 

.125 

.0008 

.0116 

.0506 

.5186 

- 3.538 

.150 

.0012 

.0127 

.0387 

.4332 

- 3.283 

.175 

.0015 

.0135 

.0289 

.3549 

- 2.975 

.200 

.0018 

.0142 

.0209 

.2847 

- 2.628 

, .225 

.0022 

.0146 

.0146 

.2236 

- 2.261 

.250 

.0026 

.0149 

.0096 

.1717 

- 1.893 

.275 

.0029 

.0151 

.0059 

.1288 

- 1.541 

.300 

.0033 

.0152 

.0032 

.0944 

- 1.220 

.325 

.0037 

.0153 

.0012 

.0675 

-.9381 

.350 

.0041 

.0153 

-.0003 

.0471 

-.7012 

.375 

.0045 

.0152 

-.0012 

.0321 

-.5093 

.400 

.0048 

.0152 

-.0019 

.0213 

-.3596 

.425 

.0052 

.0152 

-.0023 

.0138 

-.2467 

.450 

.0056 

.0151 

-.0026 

.0087 

-.1645 

.475 

.0060 

.0150 

-.0027 

.0053 

-.1066 

.500 

.0063 

.0150 

-.0028 

.0032 

-.0672 

.525 

.0067 

.0149 

-.0029 

.0019 

-.0412 

.550 

.0071 

.0148 

-.0029 

.0011 

-.0245 

.575 

.0075 

.0147 

-.0029 

.0006 

-.0142 

.600 

.0078 

.0147 

-.0029 

.0003 

-.0080 

.625 

.0082 

.0146 

-.0029 

.0002 

-.0044 

.800 

.0107 

.0141 

-.0028 

.0000 

.0000 

1.000 

.0135 

.0136 

-.0027 

.0000 

.0000 

1.40 

,0187 

.0125 

-.0025 

.0000 

.0000 

1.80 

.0255 

.0115 

-.0023 

.0000 

.0000 

2.20 

.0279 

.0106 

-.0022 

.0000 

.0000 

2.60 

.0320 

.0098 

-.0020 

.0000 

.0000 

3.0 

.0358 

.0090 

-.0019 

.0000 

.0000 

3.6 

.0409 

.0080 

-.0017 

.0000 

.0000 

4.2 

.0454 

.0070 

-.0015 

.0000 

.0000 

5.0 

.0505 

.0060 

-.0012 

.0000 

.0000 

5.8 

.0549 

.0050 

-.0011 

.0000 

.0000 

7.0 

.0603 

.0039 

-.0008 

.0000 

.0000 

8.0 

.0638 

.0032 

-.0007 

.0000 

.0000 

10.0 

,0691 

.0022 

-.0004 

.0000 

.0000 

12.0 

.0727 

.0015 

-.0003 

.0000 

.0000 

14.0 

.0752 

, .0011 

-.0002 

.0000 

.0000 

16.0 

.0771 

.0008 

-.0001 

.0000 

.0000 

18.0 

.0786 

.0007 

.0000 

.0000 

.0000 

20.0 

.0798 

.0006 

.0000 

.0000 

.0000 

22.0 

.0809 

.0005 

.0000 

.0000 

.0000 

23.6 

.0816 

.0005 

.0000 

.0000 

.0000 





















Figure 2. - Dimensionless temperature distributions for 
various Erandtl numbers. 
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Figure 3. - Comparison of small plate experimental and theoretical 
velocity distributions for Prandtl number of 0.72. 
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Figure 4. - Comparison pf small plate experimental and theoretical 
temperature distributions for Prandtl number of 0.72. 
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Figure 5. - Comparison of large plate experimental and theoretical 
velocity distributions for Prandtl number of 0.72. 
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temperature distributions for -Erandtl number of 0.72. 







Figure 0, - Dimensionless heat- transfer parameter as function of Frandtl number. 
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